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ABSTRACT 


The numerical frontogenesis model of Williams [1974] is 
Semverted to nm p vertical spacing and modified to include 
a higher upper boundary, tropopause and stratosphere, mois- 
ture and variable vertical turbulent diffusion coefficients 
of momentum and heat. The numerical solutions show more 
realistic Guasi-Steady fronts forming within 1-2 days. 
These solutions are examined and compared over a range of 
various boundary value coefficients and various ranges of 


meomperature. 
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I. INTRODUCTION 


Hoskins and Bretherton [1972] have shown analytically, 
and Williams [1967, 1972] has shown numerically that dis- 
Continuous fronts can form within a finite period of time 
if no turbulent diffusion is present. These studies sug- 
Best that a discontinuity will form within 24 to 36 hours 
when reasonable initial conditions are used. If turbulent 
Garsion 1S present, 1t can be expected that a balance will 
be achieved between the frontogenetic advections and turbu- 
Jent diffusions of momentum, heat, and moisture. The front 
should remain in this state of quasi-balance as long as the 
iamee scale detormation fleld causes frontogenetic advec- 
tions around the front. 

The numerical frontogenesis model of Williams [1974] 
etons the effects of horizontal and vertical turbulent dif- 
fusions of momentum and heat in the formation of these 
quasi-steady fronts. This model which includes an Ekman 
boundary layer will hereafter be referred to as W74 and the 
emcee aseit existed an Williams [1972] will be referred to 
as W72. Cornelius [1974] added moisture with its subsequent 
condensation and release of latent Heer GO. w 745, = iat sen 
Vestigation raises the upper boundary, includes a tropopause 
and stratosphere and appends variable vertical turbulent 
momentum and heat diffusion coefficients to the work of 


wormmelius [1974]. 


IZ 





The purpose of this study is to obtain and examine frontal 
POkmearonsmenrOUgtOUL aetOUr-day period. The physical model 
is essentially the same as one of the models treated by 
Hoskins [1971] and Hoskins and Bretherton [1972]. The 
basic model, as it now exists, numerically describes frento- 
genesis which is forced by a nondivergent horizontal wind 
field which contains stretching deformation. This deforma- 
tion wind field is constant in time and independent of height, 
except in the Ekman boundary layer. The moist hydrostatic 
Primitive equations with diffusion are used in which the 
mene ede DpenGgent quantities are functions of y and p only. 
iilemctore. sOuLpuc £10m the model is in the form of vertical 
@mess-Sections in the Y-P plane. 

ihesmodel of W74 has been transformed to pressure coordi- 
Mares £0 Obtain more accurate values of pressure for the 
Homsture equations of Cornelius [1974]. 

iimocertion 11, the basic forecast equations are developed 
and a simplification which keeps the problem two-dimensional 
ima1seussed, the initial conditions are presented in Sec- 
meoue lil, Convective adjustment processes and latent heat 
felkease in the model are included in Section IV. The experi- 
ments Undertaken are examined in Section V and conclusions 


mee efivVen in Section V1. 
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If. BASIC EQUATIONS AND BOUNDARY CONDITIONS 


The hydrostatic primitive equations with diffusion and 


moisture may be written as follows: 





dV d(wV) 
pe CUD eters Ve (RXV) = AV’ V 
, al 
bs op CK on? ’ Cen) 
Pome wc (so =] fees + 2 (x. 22 5 
at ~/ “ 3p 3 dp ‘6 Op 
Ohm Cecelia (2.2) 
3 3 9 3 
ee 8 tqv) * = (uq) = ALV'q + 5 (K, 5p) 
a, = Mowe (2.3) 
JW 
= 0) 2.4 
ae (2.4) 
ap 
ee LP) (2.5) 
where 
L(p) = = (2-)S eco ee (Rh). 
jo ee jy) ei 


iiemequations are in pressure coordinates with slight modi- 
ieeation due to the addition of the equation of conservation 
of water vapor (2.3) which is analogous to the thermodynamic 
eemation (2.2). Sinks and sources of heat and moisture have 
Been added ¢o the thermodynamic equation (2.2) and the con- 


Semvarron OL water Vapor cCquation (2.3). 
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emia tilt Les, Ayy> Ag and - ote Nee nOn Zonta turpu- 


lent diffusion coefficients of momentum, heat and moisture, 


mespectively, while Kay Ky and Sa are the corresponding ver- 
medal COCErIGlents. These coetficients were assigned various 
reasonable values in the following experiments. Dry and 


merst COnNVeccive adjustment processes are represented by 
ime Lunctions Qa and Quy while latent heating due to large 
Seate condensation is represented by Ho. faves folblenoie rian Ma 
Mena Sink Or source of moisture due to moist convective 
adjustment in the moisture equation, while Mo is a sink for 
moisture due to large scale condensation. 

Boundary conditions for the above equations are such 
that the domain has a flat lower boundary and a constant 
meessure SUrtace where w = 0 as the upper boundary. The no 
slip condition is used at the lower boundary and the no 
stress condition is used at the upper boundary. The condi- 


mens are as follows: 


w(x,y,0,t) = o(%,y,p,,t) = 0 


aie, y 5051) sent 


dV 


ogee (CSSD 5) Saw 
M op (%?¥>P,>t) (2.6) 


00 


0 ; _ 
a) op (x,y.P,,t) =a.0 


Ow 
ra) 
a 

@D 


Gey ne) = K 


K at (x,y,p,,t) yo 


p 
K es Caviar eee 
ee os? q op 


q a 
where Pp, represents the pressure level at the top of the 


domain, and 0 corresponds to Py, at the bottom. 


iS 





icmp OrsiiteematCadyactate SOlutions to the hydrostatic 


imi ervie equattens [(2,.1) 5- {2.5)] are: 
V=U =D [x(1-e °cos et ye ° Saye 
- y(1-e °cos oC) - xe ° sin i ae t2a7 ) 
WwW = Q, 
$= @ = -D*(x?+y?)/2 - £Dxy + F(p), 
uO ( Di, 
q = EN: 
where 
fe 1, i 
c= ( 7K, ) (Po°P)> Ky is evaluated at Dae and £ are 


constants and where F(p) is such that the hydrostatic equa- 
mone 2,5) 1S satisfied. The functions G(p) and H(p) are 
@aemirary Lunctions of p since ®™ = 0. These solutions sat- 
meeeeeenc boundary conditions (2.6) and allow for inclusion 
oememe Ekman boundary layer. The relations (2.7) satisfy 
mitemoccady-State Equations if some small advection terms in 
Mmmemooundary layer are neglected as they usually are in 
Ekman theory. 

The equations and the model become two-dimensional as in 
W74 and Cornelius [1974] after a subdivision of dependent 


Wemraples as follows: 


peewCo sp) + Ulty.p,t)a + vVO.p.t)) 

wo = wly,p,t) 

6 = 6(y,p,t) (2.8) 
eae ae) (yp) 


CeeedCy 5st) 
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It is noted that all departures from the approximate 
steady deformation solutions (2.7) are assumed to be inde- 
feaichteOil—~r lf we Substitute the relations (2.8) in the 


Seomponent of (Z.1) we obtain: 





du ) 3 ~ J 
are By Re or 3p (co) uD(1-e°%cos tc) - Diy(l-e S cos G) 
- xe > sin C | ha + yDe~> sin Ge a0) a [xe 5(1-cos ow) 
“¢  & _ d*u 3 du 
+ ye Sime — eye + Ay ay? + yp (Ky, 3p Ne (2.9) 


Pieeiesestiayveweswilie neglect the terms in (2.9) where 
Beappears and apply this and the other equations at x = Q. 
fiemthe initial u;, v, w, 9, n and q fields are independent of 
Peeve error from this approximation will grow slowly and 
feeeiebe Confined to the boundary layer since the x dependent 
@emns Of (2.9) are zero outside the boundary layer. Since 
magneorgenesis occurs very rapidly, it is expected that the 
development of x-variations in the dependent variables would 
mave: Only a small effect on the resulting quasi-steady front. 

iMimanywecasce these effects could not be observed in the 
menespnere, because variations in the basic fields along 
atmospheric fronts would be much more important. 


When equation (2.9) is evaluated at x = 0 it becomes 





at + (uv) + ss (GCSES) ee 7 ry ha ais aM (eS SG.) 
_ y Se Stee Bipill ego 
V yD (e Si Giles: iy + Ay Bye * 3p (Ky, “ PREZ 10) 
where V = -Dy andT = 1-e° cos erence S) component pOLe (2.1) 


Tobtedeatex = O takes the form 


a, 





OV re) 


‘) .) 
ve ay fevev in t aye Soe oer ay OY) =i By 
: ree ee ie 
(wV) ap Ce GOS) dy fu eh Ay 
dV 
(Ky ap ? 


atu, 2 
dy f) 


(asia 


Pamitlariy the following equation for the departure of poten- 


Mmeetemmerature can be derived from (2.2): 


e) Oo - 
(v6 ios 3p Coe et vy oy = Ag yy 2 


+ Qa + Qu + He 


ae, a 


E oy 





aes the following equation for q from (2.3) 


dt dy y Gey 


Zz 
24+ 2. (vq) + 5 (wq) + TV 2 = a Sos 


aa) eer I 


A 


Pemacions (2.4) and (2.5) become 


av, OL 


dy * op 
and 

oii 

op Ga) 4. 


(- 


00 


a (Ky a 


(Qa, 


es 
p M op 


(as 


(25145) 


(2245) 


ine boundary conditions at the bottom and top of the domain, 


imespeCtively, are 


= = = K CAS) = oq = = 
u Vv Ww) Ko D Ky yp 0 p Be 
. du = ov = = 00 = oq = 
eee 28 than | <q ap 


(272.16) 


~ Prop 


Piewoeecerame thesvertical average of a quantity as 
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Pp 


O 
c( ) oe — 4 ( )(ap), (ay) 
Po Ptop if 

top 


integrate the hydrostatic equation (2.15) with respect to p 


and remove the vertical mean, it follows that 


p Da (2s) 
» J a i L(p)edp - < f L(p)9dp >]. 
Prop Ptop 


fone the vertical average of (2.14) and using the boundary 


wee 1G) gives 


os 
dy 


<v> = Q. (2) 
iimeseequation states that the total mass flux in the y-direc- 
Ewen 15 independent of y. In W/2 a symmetry argument was 
used to show that this flux must be zero. This argument 
Geeemnot hold strictly in this experiment, but it can be 


Pemected that at large distances from the frontal zone, the 


meoturbance mass flux will vanish. Thus it is assumed that 
Save = 0% (Z2.20)) 


homatwon (2.11) can be written 


a + sy Ly - <vv>]j + ap) ap py (TvV-<PveV) 
; Sy [ue sin ¢ - <ue ° sin g> - Vie $5 (e°° cos £) 
- <u s (e°° cos t)>] = 7 Gis) a a 
* 5 Cy = : = ( 3 Jp=p a ay 


iitemequations (2010), (2.12), (2.13), (2.14), (2.18) and 


Cr Wjeteuiiea complete set which can be solved by a pure 


ibe 





Piieraneg process. Ihe finite difference equations conserve 
mean squares in the advection terms and are described by 
Williams [1967]. The Matsuno scheme is used everywhere 
emecpt in the arrangement of variables and the finite dif- 
fenvence approximations are the same as those used by Williams 
Meo). In order to close the problem, computational boun- 
daries must be introduced in y. Since the disturbance 
velocities should die out at a sufficient distance from the 


ms Of dilation, then 


ype 0. (2222) 


However, there is appreciable inflow across these computa- 
miomal boundaries since V{tY) = + DY. The quantities u, @ 
and q which are advected across the boundaries must be spe- 
Cified independent of the interior values if computational 
meee lity 15 to be Maintained (Platzman, 1954). Thus, the 


femehowing boundary conditions are used: 


Ue ee ay,/ 2) 3p 04) 
2) ee eer) 7 jie) yy 


Paesen + fy /2), p,t] 
per + Ay/2);, p,t] 


(2225) 
eyes y/2)5)p,t] = 8[t(Y - Ay/2),p,0] 
ql#(¥ + Ay/2), p,t] = q[+(¥ + Ay/2),p,0]. 

The computational boundaries y = +Y are placed between grid- 


ments sO that the above conditions are actually applied at 
pa <(Y + Ay/2). 

The boundary conditions were found to be satisfactory 
Peeecpt that a solution Sséparation developed near the boun- 
daries. This was controlled by using a Matsuno time dif- 


ferencing scheme for stability and because of its damping 
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@iieaeteristics. Ihe vertical gridpoints were spaced pro- 
portional to &n p using the relation p = Pau e SK 


where s is Ap/H. Refer to Table II for forecast levels. 
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TIE. INITIAL CONDITIONS 


The initial conditions of Cornelius [1974] were modified 
by raising the upper boundary to approximately the 100-mb 
Sao, and by simulating a "tropopause" centered at about 
mre 250-mb level. The Cornelius [1974] initial conditions 
mmemaeveloped as follows. The initial potential temperature 
field in pressure coordinates is given by 

_ ee p yk 
o(ysP.0) = By (py) exp [Faye - CR) 


- a(2/T)arctan(sinh ay) S209) 


where 


au a es! 1165152) 


HVL (Py) (-88,/ 3p) HY(R/Pyg) (D/P yy) * (2 07/ 2P) p 
ieee mode], a 1S given a constant median value by solving 
(ee) at p = Py amcduusine the static Stability 90, /9p at 
meee level . 

lve sinitial uw component of the velocity is found using 
peomm ot the hydrostatic equation (2.26) and the definition 


of geostrophic wind 


Et 2o 
Peesorm Of the thermal wind equation is found to be 
Eee te URS Pp 4s 20 
op feo y 10 te ap ( Py ) oy ° (3.4) 
We again assume a median pressure level Pay where UCP) = 0, 


and integrate 


OZ 





Pou Roe fl. p yk 
J pales © oy “i = ( D. )” dp (5.5) 
M M 
Bom fet 
a(y,p,0) = 2 22 po Re yk - (MK (3.6) 
yes fe 8 : Po Po . 
where 
= R : 
K [C., 


meom cquation (4.3), in Williams [1972], it is seen that 


QO 


lie eke! 
5 ee Sech (oy). (Coe) 


itemetore, the initial u field is given by 


Cc Je 
u(y,p,0) = - -p “8% (¢ B® - ( * )"}sech (ay). (3.8) 
O O 


The initial divergent wind is determined as in Williams 


wey? | by using the quasi-geostrophic equation in pressure 


Goorcinates. The stream function is defined as 
v = oy/dap 

and (3.39) 
w = - dop/doy 


and equation (4.10) in Williams [1972] takes the form 


30 
a7 | I 92 _ 2D 36 
ap2 + L (Py) C- pp, ae Fz L (Py) ay” (Sebo) 





recalling that 


Pp 
(3) = = ( a i 
iiesequatton (5.10) 1S solved for ~ from the initial temper- 
ature field (3.1) with the technique of Ogura and Charney 
[1962] and ~ is then determined from (3.9). 
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These initial conditions do not satisfy the boundary 
conditions on u, v and 6 at the upper and lower boundaries. 
Hence a period of adjustment will be required to form the 
surface friction layer and some oscillations may be observed 
later. The fields of w and m7 are obtained through the con- 
mimmuaity and hydrostatic equations (2.14) and (2.15). 

The initial field of specific humidity is calculated 
by first determining saturation vapor pressure e. at each 
gridpoint through use of an integrated form of the Clausius- 
Clapeyron equation. According to Hess [1959] 


L. 1 
lee so LI exp ( fee ( = 
V O 


| 


)) ee he 


where be Toetuiemeacenn eat Of Condensation , RY 1S the eas 
constant for water vapor, and Ts ESa yCOrstant of iIntecna- 
tion. With knowledge of the saturation vapor pressure, 
Mremsattiration Specific humidity is obtained from the ap- 


weoecimate relation of Haltiner and Martin [1957] 


BOLZ e. 
| as eae (3212) 
p 
Knowledge of potential temperature at t = 0 from equation 


fot), makes it possible to calculate q. at each point using 
cametrons (3.11) and (3.12). Defining q to be a given 
ieaction of qd. ablewsethe testing of atmospheres of different 
womstant initial relative humidity. The initial conditions 
of Cornelius [1974] were modified as follows. A tropopause 
mas Simulated centered at about the 250-mb level. This was 
Aeconiplrshed Wsime a reference temperature at the midpoiit 


of the horizontal scale and calculating temperatures such 
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that a tropopause with its expected variation in height was 
feeecrmined and a Zero vertical temperature gradient was 
achieved in the lower stratosphere. The tropopause height 
[err caeerom o25 to LlQ mb over the horizontal domain for 

8, @,) = 290 K. The tropopause was simulated as follows. 

Mm reference temperature was computed using equation (3.1) 
as 


js © 
It = 215 +a ( & ) arctan(sinh ay) ( pO (Scho) 


ia 
Ee ins 
hirem the temperature in each column reached this reference 
temperature, it was held constant at that value, as expected 


in the lower stratosphere. Temperature was then converted 


back to potential temperature. 
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im, CONVECTIVE ADJUSTMENT, LATENT HEAT AND PRECIPITATION 


Two additional heating terms Os and He) were included 
tmeene moist thermodynamic equation (2.2) necessitating 
further amendment to W74. Atmospheric heating by the dry 
convective adjustment, Q? was included in W74 to eliminate 
Giyeititense grid scale convection which developed due to the 
Mmemiation Of any dry unstable lapse rates. This redistribu- 
tion of heating was accomplished using averaging techniques 
to ensure that the dry adiabatic lapse rate was not exceeded 
fecyop < 0). 

Lapse rates were tested in each vertical column of the 
model as in Cornelius [1974]. Upon encountering an unstable 
layer (d30/dp > 0), a cumulative average was taken and applied 
Memmme Sridpoints within this layer until neutral stability 
(90/dp = 0) was attained. The model was checked simultaneously 
for instabilities developing below the convectively adjusted 
Mayer and potential temperatures at these levels were also 
cumulatively averaged and reassigned the new values. The 
final result of the technique is that the heating is redis- 
tributed to simulate physical processes of the atmosphere 
and mathematically that 00/dp< 0 throughout each column of 
the model. A completely analogous scheme was used for satu- 
rated parcels in the moist experiment, only using equivalent 
merential temperatures, Bop» Inbaece Of potential temperature 


0. 
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iiewidd@tton or moisture required that large scale con- 
densation and the subsequent release of latent heat be given 
miportant consideration. Following Haltiner [1971], changes 
mmemersture and temperature due to large scale condensation 
were determined as follows. 

Where relative humidity exceeds 100%, the adjustments in 
temperature 6T and specific humidity 6éq are determined by 


ene following Equations: 


cme) 7 Olt; )) (4.1) 
and 


c,éT = - L.éq. (Ae) 


iitese Cuuations require that the final state of the air be 
exactly saturated and that excess moisture be condensed iso- 
iioically, releasing the latent heat to the air. 

Hauatsons, (4.1) and (4.2) can be expressed in iterative 


temm as 


aq 
eee *) 1 7 ), ANd = ra. (4.3) 


and 


CAT = Eu Ne) (4.4) 


fminating AT from (4.3) and (4.4) gives 


oO 

ea TE ey es 
1h 
4 r) 


At constant pressure the Clausius-Clapeyron equation can be 


expressed as 


q 
=r Jy = RT TE: (4.6) 








Pumstituting (4-6) into (4.5) gives 


q,(1 - 1) ; 
Aq = = ei a (4.7) 


P+ 2 
Cah! 





Aq may be now calculated from (4.7) as can the corresponding 
Meeuic of AT from (4.4). The first estimates of the adjusted 


Bemperature and specific humidity are, 


J I eel and Ge ge AG (4.8) 


These values in turn are used to obtain improved approxima- 
trons. Iterations are continued until the values of T' and 
Geet Cach level match with sufficient accuracy 100% rela- 
mievcenumidity. dhe changes are then applied to the potential 
Mamiperature and specific humidity fields of the model ac- 
Counting for the term H of equation (2.2) and M. of equation 
5). 

The model is then investigated for unstable moist convec- 
megnmand itS associated heat of condensation. A convective 
adjustment scheme similar to the dry one is used on the satu- 
rated parcels so that moist static stability is always main- 
tained. A field of equivalent potential temperature is 


Bomnibated according to Holton {1972] as 





Hews 
Oc1; = 0 exp — 7 ) (4.9) 


and the moist convective adjustment is designed to ensure 


mat the moist adiabatic lapse rate is not exceeded 


(30.),/9p Se 
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iitewicat Trehease . oe Poa result Or MOrst Convective 


adjustment is included through knowledge of the change of 
pauivalent potential temperature during the adjustment. 

Taking the natural logarithm of both sides of equation 
(4.9) gives 


Liq 
CO meaner CO a 


SE ommgee bi (4.10) 
p 





A constant pressure surface may be assumed since the model 
Mees pressure coordinates. Differentiating both sides of 
(4.10) and recalling the definition of potential temperature 
gives 

en _ aut 


ae aa s — + LTT Se | aaa i» Caramamar e (4.11) 
one Tr a i T 


Solving equation (4.6) for dq. SUbDSELLULING Into equation 


(4.11) and solving for dT leads to 


SOE 
ce ae ye ee els) 
Se. Rec) 4 cual 
v Pp Pp 


where dOop PSethewenamge in Ocr needed to adjust (90o,/9p = 
De 

Since the atmosphere remains saturated throughout the 
MeeGess Of moist convective adjustment, the Clausius-Clapeyron 
equation is again used to obtain the corresponding dq or 
actually dq... Chamees ware sagaim applied to the potential 
memperature and specific humidity fields of the model account- 
ing for the terms Qa Gt eaGUctEtOtew2. 2) and ne of equation 


2. 5) . 


Zo 





A cumulative amount of precipitation is then found using 
the relation 


Po 


: éqdp 


O 


GQ |b: 


(Axis) 


mreagreement With the changes in specific humidity due to 


condensation and moist convective adjustment. 


Precimitae won 
ieee xpressed in inches of rainfall in the model. 





V. NUMERICAL SOLUTIONS 


All numerical results to be shown use the following 


wemues for the constants: 


meo= S40 sec Py = ONO sede Puy = See 
wey = .01978 Y = 1800 KM Ay = 60 KM 
g = 9.81 msec f= 10° sec D= 10° sec 
2 = 300 K 8, (p,) =a .2 20,500) K Cy = ae ee 
Eee - 461 joules R = 287 joules Kg K™’ 
ih coal 
Mg K = R/C, L. = 2.5x10° joules 
, 7] 
a = 12.56 K RH = 70% Kg 
08 08 
I = r =e I Eon 4 a -1 
a 5S K km (“ap py =—-55446"K cb 
A, = A, = 3x10*m’sec * A, = 100000 m?sec — 
rae | = 
€, = 0, cb*sec Cap = 0, .001409, .0007 cb*sec — 
(which corresponds to 0, 10, S m?sec. in hevont coordinates) 
Con soy 20007, . 0001761 cbh2sec (which corresponds to 


0, 5, 1.25 m¢sec.” in height coordinates) 
where more than one value means that the constant was varied 
Mmoimrerent experiments. Whenever possible, the values are 
the same as those used in W74. 

imtes Values of vertical diffusion of moisture Cs and heat 
Cy Meme Varled in the experiments to compensate for density 
fewer lon 1n the atmosphere. This is equivalent to a constant 
meeteiretent in the z coordinate system and is accomplished 


using the following formulac: 


oa 





Coal (ile Pee (5.1) 
and 

coche 0/ Da) s C553) 
in this section numerical frontogenesis solutions with mois- 
miiee, friction and a stratosphere are examined at t = 4 days 
and compared to the corresponding dry atmosphere. 

ro cam De Seen trom Table 1, 11 experiments were performed. 
meeriments 1 and 2 contain only horizontal turbulent diffu- 
Sions of momentum, heat and moisture while Experiments 3 
through 11 contain both horizontal and vertical turbulent 
diffusions of these quantities. 

In Table I, 8, (p,) is the reference potential temperature 
(K). The condensation interval is 5 time steps between suc- 
cessive applications of condensation and moist convective 
adjustment to the model in the moist experiments. Horizon- 
Pemedrbulent diffusion coefficients and initial boundary 
values of vertical turbulent diffusion coefficients (Cp: 
Cop: C.? are relative humidity (percent) and initial static 
Seaplity (K/cb). 

Peaperinenits were performed with a reference potential 
temperature (0; (pP.))> OLezsuky 2O0KeanG S00K to simulate 
pieereally realistic situations. Changes to the dry atmo- 
mmeere With different reference temperatures were almost 
feeieecible, but due to the exponential function in the 
memsius-Clapeyron equation important differences were ob- 


Served in the moist atmosphere. 


oy 








Eaepelimemts were also conducted using different boundary 
values of vertical turbulent diffusion coefficients of momentum 
and heat. 

The final experiment included a turbulent mixing coeffi- 
Cient which reduced the values of Cap and Cop linearivaieo il 
= their original values between approximately the 93- and 
ef-cb levels in five equal steps. 

ine adjusted values "of vertical turbulent diffusion co- 
erracients of momentum Cs and heat Cy, were then varied as 
@eceribed by equations (5.1) and (5.2). 

Primary emphasis on analyses of experimental results was 
Gieecontour plots from the IBM 360/67 computer. These results 
obtained for all 11 experiments are analyses at t = 4 days 
G@eerme fields of u, q, 8 and w. The entire vertical plane 
meeoown irom 100 to 10 cb in these plots. However, the 
outer portion of the domain between |y| = 1200 km and ly] = 
rovekm 1s not shown. 

Peeures 1-3 show the initial fields of u, 8 and w at 
meeerence potential temperature of 290 K used in the follow- 
micnexperiments. Contour plots for experiments 1 and 2 are 
shown as figures 4-10. Figures 4 and 5 show distributions 
of the rotational wind, u for corresponding dry and moist 
experiments, respectively. Features of u at all levels dif- 
memeeguite significantly from dry to moist experiments. There 
is a marked increase in horizontal and vertical shear in the 
easterly flow in the moist experiment in the vicinity of the 
Hront near the surface. Horizontal and vertical wind shear 


momalse more intense at low middle to middle levels in the 
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Mest ecxperiment. If is noteworthy that this property is 
dominant throughout the moist solutions. In the moist ex- 
periment the jet core increases in magnitude while shear 
decreases to zero at a lower level in the northern region 
Mear the tropopause. Figures 6 and 7 contrast values of 
6 throughout the domain. Noted is the warmer air in middle 
ewels in the frontal zone. A slightly more intense gra- 
dent iS seen near the surface with little change in tilt 
Gmmeme front. Ali this is attributed to the latent heat 
memsease in the moist experiment 2Z. 

wertical MOtlONs fOr the experiments are investigated in 
Figures 8 and 9. Upward vertical motions occupy a much 
Pageer zone, extend much higher and occur closer to the 
meomcal zone in the warm air mass of the moist experiment. 
Significant downward motion occurs between |Y| = 1200 and 
maie= 1800 km (not shown in the figures). 

meine lO) sShews the field of specific humidity for ex- 
Pommmemt number 2. Maximum precipitation is .37/ inches. 
Pieenes Il through 17 compare the dry and moist fields of 
Uso) and w and show the q fields of experiments 3 and 4. 
Figures 11 and 12 show a small increase in the zone of 
feee@emty Component in the lower levels of the moist experi- 
fee vertical and horizontal shear are also seen to in- 
G@mease aS in experiments 1 and 2. Also noted is a much 
larger jet core in the moist experiment, with a slight 
iemunwerd shift. Noteworthy is a significant decrease in 
ieeeitude Of the casterly component of u near the surface 


Pome tiat Of experiments 1 and Z. Figures 135 and 14 exhibit 
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tiieesame COntrast in potential temperature as that of ex- 
meamments 1 and Z. If 1s noted that the air at low levels 
ficenear the tropopause 1s colder in the moist experiment 
fmen vertical turbulent daffusion added (exp. 4). This is 
ase truce in the dry experiments. This is attributed to 

the increased amount of heat loss due to the vertical tur- 
bulent diffusion of momentum and heat offsetting the release 
Gelatent heat 1n the moist experiments. The air at mid 
Tevels in the frontal zone 1S warmer in the moist experi- 
ment as seen before. Figures 15 and 16 show a marked change 
immevertical motion between dry and moist experiments. A 
Harpe zone of upward motion is seen near the frontal zone 
With maximum values exceeding 3x10 mb/sec near the 85-cb 
Memel. Figure 17 shows the specific humidity of experiment 
MEO LeWOrthny 1S the dramatic shift to the south in specific 
homadtty following the frontal position. 

Faetures 19 through 38 compare the moist and dry fields 
Seo and w and show the specific humidity of the moist 
EememamentSs using boundary values of vertical turbulent dif- 
fusion of heat and momentum of .0007 cb2sec. (5 m’sec’’ 
memieront coordinates) at reference potential temperatures 
eieeeO, 280 and 300 K. 

mertnes 1s and 19 contrast the w tield of experiments 5 
and 6 at 8, (,) = 290 K. Noteworthy is the increase in zone 
of maximum easterly component of the wind at low levels and 
the more intense horizontal and vertical shear in low levels 
Pyeoreexperiments 3 and 4. A slight northward shift and in- 


Smease 1M size Of the jet coré™is noted in the moist experiment. 
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Peepunes, 20 and 21 show little difference in potential tenm- 
Wemaeure ficlds, although the air in the frontal zone in 
middle levels is warmer in the moist experiment due latent 
Meremrclease. Ihere is a marked increase in tilt of the 
Meee OVEr previous experiments, due to more normal north- 
Mera position of the front at the surface. Figures 22 and 
Momeomcrast the vertical motion of experiments 5 and 6 show- 
ing a maximum in upward motion near the 850-mb level con- 
Sweet With previous experiments. Figure 24 shows the 
Seeemric humidity field which differs slightly more in the 
Paiiesmanner as experiment 4 differed from 2. Maximum pre- 
@mpictation is .4 in. 

Pesirese zo and 26 show more intense horizontal and ver- 
meerleshear apparent in the easterly component of the moist 
experiment at 8,(p,) = 280 K. The jet core is reduced in 
eeemand magnitude from experiments 5 and 6 and the hori- 
momma shear decreases significantly near the tropopause. 
neeeamoulad be noted here that the tropopause in experiments 
7 and 8 is lower than the tropopause in previous and subse- 
Weememexperiments due to the lower reference potential tem- 
Pemeeire used. Figures 27 and 28 show little difference 
wemmecn moist and dry fields of potential temperature except 
femeeene warmer air in mid levels in the frontal zone of the 
iemste EXperiment due latent heat release. Noteworthy in 
these experiments 7 and 8 is the zone of low stability and 
meperent break in the tropopause in the northern part of the 


region. Figures 29 and 30 show a similar contrast between 
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itomsG and dry fields of vertical motion as seen in previous 
eepetimentsSs except that the contrast is not as dramatic as 
seen before. Figure 31 shows the q field is reduced con- 
siderably over previous experiments due to the lower refer- 
ence potential temperature used (280 K). Maximum precipita- 
mionmereduced to .15 inches. 

Figures 32 through 40 show a much higher tropopause 
aiemeo the larger reference potential temperature (300 K). 
A dramatic change in the u field is seen in Figures 32 and 
33 between dry and moist fields. Especially noteworthy 
7S the very strong jet core and much more intense horizon- 
tal and vertical wind shear. Also note a large shift in 
tiemyet tO the north which was not quite so obvious in pre- 
Vious experiments. Figures 34 and 35 also show a marked 
contrast between dry and moist potential temperature fields. 
The air iS much warmer in the upper levels near the tropo- 
pause due to the latent heat release, Thermal gradients 
fiomeeene £trontal zone are also stronger in the mid and up- 
Memecvels, Strong heating in upper levels 1S apparent. 
femtneal motions {Figures 36 and 37) exceed 3x10 “mb/sec 
Mmmeme Warm MO1St air. The contrast 1S similar to that seen 
mere vious experiments. The moisture field (Figure 38) is 
greatly intensified over previous experiments and specific 
humidities are in excess of 19 g/kg at low levels to the 
Seuem, This 1s no doubt due to the larger quantities of 
latent heat release observed in this experiment, and the 
Pier imtttal reference potential temperature. Maximum 


Puceipitation increased to .71 inches. 
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Figures 39 and 42 show experiment 11 which included the 
turbulent mixing coefficient into experiment 4. The only 
fieererences noted were a small pocket of lower specific 
imeatey (Freure 42) at the top of the mixing level and a 
decreased zone of maximum upward motion in the warm moist 
aap@eoree the mixed experiment (Figure 41). From the above 
Womeomelude that vertical turbulent mixing has no effect 
on mid and upper level deformation fields. 

Comparison of the results of the previous experiments 
with similar experiments of Cornelius [1974] revealed the 
following: 

1. The westerly wind field extended much lower while 
mitemedsterly wind field was compressed below the 85 cb 
ieveiwee Lhe zero wind level was also observed to be lower 
tee CXpanded atmosphere with tropopause included. fhis 
Pec in both moist and dry experiments. The jet was also 
lower and more apparent. 

Meine pockets of warm alr were apparent in the frontal 
zone as in Cornelius [1974], but frontal tilt was somewhat 
less after 4 days in the experiments conducted in this study. 

Ewe olegniticant downward vertical motion occurred outside 
the zone observed in both studies although larger values are 
evident in the studies of Cornelius [1974]. 


a specitme humidity frelids were similar in both studies. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


The numerical model of Williams [1974] is modified in 
mais Study to include moisture, a tropopause and lower stra- 
tMeopmere and variable vertical turbulent diffusion of heat 
and momentum. This leads to a more realistic situation and 
Sousequently more realistic solutions. Intensification of 
miesbDaroclinic zone 1S seen to occur in upper levels, al- 
miemeh mo effect is apparent on the surface due to the in- 
SmisiOn Of mo1sture. The effects of vertical turbulent 
meetusion 1S apparent throughout the frontal zone through 
amereased tilt, reduced heating and reduced values of wind- 
speed. The expanded boundary and inclusion of the tropo- 
pause along with a conversion to &n p vertical spacing 
meameed to Shift all features to a lower level in the domain. 
wet cores were generally larger and vertical motions in- 
@meased in the warm moist air near the frontal zone due to 
the heating. The moist models illustrate an exponential 
thermal dependence in the various experiments at the reference 
Memential temperatures of 280, 290 and 300 K. 

Mimedtune studies using this model, attention should be 
directed toward the use of still more realistic boundaries, 
eulem as a higher upper boundary. A translational motion 
should be added to the system so that warm and cold frontal 
eeiectrs could be investigated. Finally, an improved convec- 
me scheme such as that of Arakawa and Schubert [1973] 


Eomlds be applied to the model. 
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Table II. Forecast Levels (cb) 


EVEL p (ch) p(cb) Diced) 
1 95.86 19 44.81 Sy 20.95 
2 91.9 20 42.96 38 ZOOS 
5 Sone zal 41.18 39 1925 
4 84.45 ae. 39.48 40 Leoe2S 
5 80.96 ZS 57280 41 7 Oo 
6 (7 ol 24 50226 42 foes 
7 74.4 ZS 34.78 43 Toe 26 
8 US S:Z 26 535. 54 44 52959 
3 GS5.57 Z7 31.96 45 14.94 
10 O54 28 30.64 46 1452 
11 G28. 5 29 js Se Si 47 ils ay A 
eZ G0n 25 30 DS ALS 48 T3516 
13 Sod eal 26.99 49 1262 
14 Ho eeiters or 7A hts) 50 20 
15 53.06 5g 24.81 51 Rs o.0 
16 owes / 34 7. yes a a he 
ey 48.76 55 22.80 Sus. LOS 6 
18 46.75 36 71, 85 54 Qe 2 
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